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Abstract 

The search for regular black hole solutions in classical gravity leads us to consider a core of 
Euclidean signature in the interior of a black hole. Solutions of Lorentzian and Euclidean general 
relativity match in such a way that energy densities and pressures of an isotropic perfect fluid form 
are everywhere finite and continuous. Although the weak energy condition cannot be satisfied 
for these solutions in general relativity, it can be when higher derivative terms are added. A 
numerical study shows how the transition becomes smoother in theories with more derivatives. 
As an alternative to the Euclidean core, we also discuss a closely related time dependent orbifold 
construction with a smooth space-like boundary inside the horizon. 
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I. INTRODUCTION 



The singularity problem of general relativity (GR), as embodied by the Schwarzschild 
(Schd) solution, may not be a general consequence of a general classical theory of gravity. 
It may merely be an artifact of truncating a derivative expansion at the lowest term, the 
Einstein action, and restricting oneself to solutions of that truncated theory. Clearly, before 
the singularity at r = is reached the Einstein action can no longer be trusted. From this 
viewpoint it is of some interest to see whether the singularity problem that arises in classical 
gravity could also be resolved within classical gravity. The form of such a resolution in 
higher derivative theories could also have a bearing on how singularity problems are to be 
addressed in quantum theories. 

We shall refer to some of our results as robust if they are likely to hold for a general 
classical theory of gravity with any number of derivatives, with the possible exception of 
isolated cases in the space of theories. These results have been found to be true for theories 
with various number of derivatives (including theories with up to six derivatives), and al- 
though we offer no proof, we expect such results to be true quite generally. Any particular 
theory involving a truncation to some number of derivatives can of course yield spurious and 
unphysical results (for example the ghosts and tachyons that can appear upon linearization 
Q]), but the set of robust results applying to the whole set of higher derivative theories 
should have more significance, and should have some bearing on questions involving strong 
gravity. 

We expect a regular (nonsingular) uncharged black hole to involve some smooth matter 
distribution in its interior. Interesting regular solutions must arise from physically acceptable 
energy- momentum tensors T^ v . In the Schd coordinate system the metric of interest is 

ds 2 = -B{r)dt 2 + A{r)dr 2 + r 2 {d6 2 + sin 2 6d<p 2 ). (1) 

In some other attempts to find regular black hole solutions in GR, a constraint on the metric, 
A(r)B(r) = 1, is imposed and T^ v is allowed to have whatever spatial anisotropy is needed to 
satisfy the equations [2j. But such assumptions are not correct for physical stellar solutions, 
and so it is questionable that they would apply to regular black holes. We are moving in the 
opposite direction in this work. In particular we will not constrain A{r) and B(r), and we 
will take T^ u to have isotropic perfect fluid form. In this framework, both in GR and higher 
derivative theories, we will see whether various physically motivated energy conditions can 
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be satisfied. Energy conditions play a central role in various no-go [3J theorems in GR, but 
see also j^J on the possible violation of energy conditions. 

If the curvature invariants R R„ P R^... are to be regular at the origin, then 

this restricts the behavior of A(r) and B{r). For small r we must have |5( 



A possible multiplicative factor in B(r) arising from a rescaling of t is suppressed. Linear and 
cubic powers are not allowed. The scalar curvature, for example, at r = is R = 6(a 2 — 62)- 
We find that A(r) and B(r) of this form provides an example of a robust solution near 
the origin. Note that by requiring a regular solution at r = we are assuming that the 
origin actually exists. Thus we do not consider the possibility that a black hole contains 
a wormhole, a horn or a flux tube [3J. But in section 4 we will discuss a smooth orbifold 
construction that yields a black hole without a center. 

In GR all the parameters in (j2J) are determined for given energy density p(r) and pressure 
p(r) profiles. This is not true for a general 2 + 4 derivative theory, with action 



In this case the coefficients a 2 and b 2 are free parameters. This freedom is required to match 
these solutions to the physical large r weak gravity solutions (assuming no singularities 
in between). This is because the linearized 2 + 4 theory has modes both increasing and 
decreasing exponentially with large r and the adjustment to zero of the two increasing 
modes requires two degrees of freedom. This ability to match the small r solutions with 
the physical large r solutions persists in even higher derivative gravity. For example in a 
general 2 + 4 + 6 derivative theory (which includes terms cubic in the curvature) solutions 
near r = exist for which four coefficients, a 2 , 6 2 , 04, 64, are free parameters. This freedom 
is required to adjust to zero the four exponentially growing modes at large r. 

We remind the reader that the Schd solution itself is not robust since it does not exist in 
general gravity theories with six or more derivatives. In p it was observed that for these 
theories, (J2J) is the only possible solution of the equations when expanded about the origin. 
This applies when A(r) and B(r) have any power law behavior near r = 0, with smooth or 
vanishing matter distribution up to a possible 5-function source at r = 0. 



A(r) = 1 + a 2 r 2 + a 4 r 4 + ... 
B(r) = ±1 + b 2 r 2 + b 4 r 4 + ... 



(2) 




(3) 
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II. A CROSSING POINT 



The important point is that independent of what the theory is, A(0) = 1 for a regular 
solution. This is a robust result. But inside the event horizon at r h , A(r) is negative, due 
to the change of signs of A(r) and B(r) at the horizon. Then to match onto a regular 
solution with A(0) = 1 another sign change is required. In other words A(r) must have an 
even number of sign changes between r = and r = oo, rather than the single sign change 
implied by a single horizon. Thus for a regular black hole with a center A(r) must change 
sign at some radius r* which we refer to as the crossing point, where r* < r h . 

We shall investigate this crossing point by performing a series expansion of the equations 
of GR and higher derivative theories. Basic properties of this crossing point turn out to 
be robust. Later we shall use numerical analysis to construct solutions for all r that are 
compatible with what we know about the origin, the crossing point, and the horizon. 

As far as the matter distribution is concerned we shall consider the physically most 
conservative form for the energy-momentum tensor in the black hole interior: the perfect 
fluid form with isotropic pressure. Given that r and t have exchanged their space- and 
time-like roles for the region r* < r < r h we define p{r) and p{r) as follows, 

T„ = pg rr , T tt = -pgtt, T ee = ~P9ee, T H = -pg^. (4) 

This isotropic form would also arise from a general scalar field action with minimal coupling 
to gravity, given that the scalar field only depends on r. But the curvature tensor itself 
distinguishes between the spatial t direction and the angular directions, and so the isotropic 
perfect fluid form could be modified through general scalar-gravity coupling terms. But 
this is not sufficient to justify the use of far from isotropic pressure in regions where the 
curvature is not large. For definiteness we discuss the case of isotropic pressure everywhere, 
in which case there are three equations (the rr, tt and 99 equations) for the four functions 
A(r), B(r), p(r) and p(r). Our basic strategy will be to specify one of these functions and 
then determine the other three, but this strategy cannot guarantee that both A(r) and B(r) 
are completely regular. Nor can it guarantee that p(r) and p(r) are physically acceptable 
and satisfy, for example, the dominant energy condition \p(r)\ < p(r), or at least the weak 
energy condition, p{r) > and p{r) +p(r) > 0. 

We have already seen that A(r) must change sign at some point, while B{r) may or may 
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not change sign at this point. Our analysis of the series expansions of the equations around 
the crossing point finds two types of regular solutions. 

Type 1) 1/A(n) = and B(n) = 

Type 2) l/^r*) = and -B(r^) is finite and nonzero 

This result is robust, and we will also find robust constraints on for each type of solution. 

The first type of solution includes the original event horizon, and so to emphasize this we 
replace r* with to discuss this case. The curvature invariants do not diverge at even 
in the presence of matter, but we find the following constraints, which thus apply to both 
the original horizon and a possible "interior horizon" . In particular 



P(r h ) = -p(r h ), (5) 

and the approach to the horizon from the two sides is constrained as follows, 1 

p' ( r h) = — 3p' (rh) when t is time-like, (6) 

3p' (?"h) = —p' (^h) when r is time-like. (7) 



When these derivatives do not vanish, the weak and dominant energy conditions can still be 
satisfied in the vicinity of the horizon as long as p{r^ > 0. For the same reason p(r) must 
increase as one moves away from the horizon on the time-like-t side, and decrease on the 
time-like-r side. We also note that if p(rh) and p(rh) and their first n — 1 derivatives vanish 
then we have the following relation among their nth derivatives. 2 

p( n ) (r h ) = -(2n + l)p (n) (r h ) when t is time-like (8) 
(2n + l)p (n) (r h ) = -p (n) (r h ) when r is time-like (9) 

These robust constraints on p{r) and p(r) are acceptable for the original horizon, where 
we can assume that p(r) and p(r) vanish, but they prove to be quite restrictive for solutions 
in the interior of a regular black hole where we expect a nonvanishing matter distribution. 
In fact we were unable to realize a solution of this type in our numerical analysis. Note 

1 On the side where t is time-like, T tt = pg tt , T rr = -p r g rr , Tee = -Pegee, = -peg<p<p- The p appearing 
in © and JSJl is actually p r . pe is completely smooth across the horizon. 

2 We don't bother to display the more general relations here. 
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deal with this see 



a, Liic 

0,3 



that if anisotropic pressure is allowed then analytical solutions with an interior horizon can 

Qfi n 
. But there are known stability problems [y], if not a no-go theorem [7[, 

associated with this type of interior horizon. 

We thus turn to the second type of regular solution where B(r*) 7^ and only one metric 

component is changing sign. In the situation that A(r) does not change sign again before 

reaching r = (if it did it would have to change sign an even number of times) then B(r) 

cannot change sign for r < r*. This is because any point where B(r) and not A(r) changes 

sign has a curvature singularity. Thus for < r < r*, g rr is once again space-like and g tt 

remains space-like, and we have a metric with Euclidean signature. This will bring up some 

issues concerning a change in signature to which we return shortly (for other attempts to 

1). 

For this type of solution we do not find robust constraints on at the crossing point. 
This is an indication that such solutions will be easier to realize numerically. But we do find 
robust constraints at the origin of the Euclidean core where 

p(0) = -p(0), (10) 
p / (0) = p'(0) = 0, (11) 
P"(0) = ~/(0). (12) 

This shows how the geometry approaches Euclidean dS or AdS space near the origin. If p(0) 
and p(0) and their first 2n — 2 derivatives vanish then we have the following relation among 
their (2n)th derivatives (assuming that the higher powers in (J2J) are even). 

P (2n) (0) = -^-rP i2n) (0) (13) 
n + 1 

We have seen that p = —p, or in other words oc g^, must hold at the original horizon, 
and now we see that it must also hold either at the crossing if it is of type 1, or at the origin 
if the crossing is of type 2. 



III. SOLUTIONS CLOSE TO THE CROSSING POINT 



We find robust regular solutions of type 2 of the following form. 

Mr) = E^-^ 



oo 

B(r) + 1 = J2 b ^ r ~ r *y ( 14 ) 

i=l 

There is some number of free parameters depending on the number of derivatives in the 
theory. If p(r) is given then p{r) and its derivatives at r = r» are determined. For GR all 
the Oj and 6j coefficients are determined, for 2 + 4 gravity a\ and b\ remain free parameters, 
for 2 + 4 + 6 derivative gravity a\, bi, ci2, 62 are free, etc. The number of parameters is 
the same as for the regular solutions at the origin in (J3J), and it is the same as the number 
needed to match onto the Newtonian solution at large r. 

The existence of these robust solutions is quite remarkable since it implies that all quan- 
tities remain finite and are smoothly matched at the boundary between spaces of different 
signature. In particular the extrinsic curvature vanishes on both sides of the boundary. We 
emphasize that we are dealing with solutions of different theories, since the \/—g that ap- 
pears in the action of the theory with Lorentzian signature must be replaced by y/g in the 
action of the theory with Euclidean signature. But after accounting for the change of sign 
of A(r), the field equations are the same in the two theories, and this allows the solutions 
of the two theories to match. 

One might wonder whether there is a single theory that describes both spaces including 
the crossing point; for example one could consider the theory with \J—g — > y\g\- But the 
theory then becomes nonanalytic in the fields, and if taken seriously, would generate various 
5-function singularities in the field equations at the crossing point. But these singularities are 
not unique, since different nonanalytic extensions of the theory can give different singularities 
while keeping the theory on either side the same. 3 This is related to the controversy in 
the attempts to define junction conditions at a signature changing boundary (see the last 
reference of js| ) . We are resigned to the fact that the theory is not well defined at the 
crossing point, and the field equations are ambiguous at this point. We will comment more 
on this at the end, but for now we will simply define p(r) and p(r) by their approach to the 
crossing point, using the theories on either side. It is significant then that solutions exist 
where all physical quantities defined in this way are finite and smooth in the vicinity of the 
crossing point. 

3 For example if e is the sign function then e(g) 2 = 1 up to the existence of singular, ill-defined derivatives 
at g = 0. 
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For these solutions (and the more general ones below) we have the following result that 
is specific to GR. 

8ttG 



r 



2 



p (r .) = ^(a.r.-l) (15) 

^ * 

1/A(r) is changing sign such that ai < 0, and therefore p(r*) + p(r#) < 0. Thus the weak 
(and null) energy condition cannot be satisfied in the vicinity of r* in GR. This is the result 
that forces us to consider higher derivative theories; higher derivatives introduce the higher 
order coefficients of (JTljl into the results for p(r*) and p(r*). It is then possible, as we 
explore below, that a more sensible T^ v can exist when the crossing point is in a region of 
high curvature. 



IV. ORBIFOLDING INSTEAD 

As an alternative to signature change, the form of the solution in (j!4|) suggests a way 
to remove the Euclidean core region and obtain a regular black hole without a center. A 
smooth orbifold around r = r* can eliminate the core in the following way. In terms of the 
proper time u from the crossing point, where u 2 = r — r*, the metric can be rewritten as 

ds 2 = -B(u 2 )dt 2 + AA(u 2 )du 2 + (u 2 + n) 2 (d6 2 + sin 2 6d<J) 2 ) 
1 1 00 

A(u 2 ) Airf + r^u 2 ^ 1 

oo 

B[u 2 ) + l = B{u 2 + n) + 1 = b i u2i - ( 16 ) 

i=l 

Recall that a\ < so that both A(0) and -B(O) are negative. The time coordinate u covers 
only r — > and the metric contains only even powers of u. We thus impose an orbifold 
condition at u = by identifying positive and negative values of u. This gives a smooth 
space-like boundary at u = where curvature invariants are clearly regular and extrinsic 
curvature vanishes. This would be an illustration of how a smooth time-dependent orbifold 
{u being time-like) resolves the space-like Schd sing ularity, see Fig. (JTJ). Related ideas have 
been receiving some attention in string theory [lCj, especially in regards to the space-like 
cosmic singularity. 
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FIG. 1: We compare the Penrose diagram for a Schd black hole (left) to the orbifold construction 
(right). In the latter, time ends on a smooth space-like boundary rather than a singularity. 

Our results concerning the weak energy condition will also apply to the orbifold construc- 
tion. There is the additional question of how a gravitational collapse of a smooth matter 
distribution can change topology to form a black hole without a center. The orbifold possi- 
bility lies outside the main development of this paper, but we will return to it briefly in the 
Discussion. 

V. MORE GENERAL SOLUTIONS 

If the solutions in (JT4"j) were the most general possible for a crossing point, then we could 
proceed to look for completely regular solutions by numerically exploring how the regular 
solutions at the origin, the crossing point, and the horizon can be joined together. But we 
find that the solutions in ()14)) are not the most general; there are more general solutions 
that are not completely regular at the crossing point. Fractional powers are involved (and e 
is the sign function). 

_L = J2 \r-r,\ l (at + a,7 e (r - r*)) 
1 } Mi,|,2,|... } 

B(r) + 1= Yl \ r - r *HK + K< r - r *)) ( 17 ) 

i={il,f,2,|...} 

This contains the terms in ()14|) since we can identify etc. 4 These 

solutions are robust, and it is these more general solutions that will typically emerge in the 

4 We implicitly replace |x|e(x) by x, etc. 
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numerical analysis. We stress that despite the fractional powers, the curvature invariants 
(and the extrinsic curvatures) continue to remain finite as the crossing point is approached. 5 
General values of the parameters in (|17|) can result in discontinuities in physical quantities. 
We will therefore study a constrained class of solutions that satisfy the following criteria. 

• Physical quantities are continuous and are as smooth as possible at the crossing point. 

• p(r) is completely regular everywhere. 

• The solutions are realizable numerically in complete black hole solutions. 

The second criteria is convenient for the numerical study, since the two field equations 
involving p(r) allow both A(r) and B(r) to be determined from p(r) alone. The third 
criteria will keep us from the completely regular solutions in (|14j). at least for the time 
being. 

We consider GR first. Here the structure of the equations imply that a completely regular 
p(r) implies a completely regular 1/A(r) of the form in ()14j). p(r) is then at least continuous 
at the crossing, and to ensure that the equations are consistent with a continuous p(r) 
requires that 

K = -UbJ 2 + bl? 2 ). (18) 



4 ^1/2 t »y 2 ^ 

This then implies that the curvature invariants R, R^R^ and R fJlU apR >J ' ue7p are also continu- 
ous. Next we must eliminate possible \r — r*] 1 / 2 or \r — r*| 1 / 2 e(r — r*) terms in the equations 
for consistency with a regular p(r). In other words p'(r) must be continuous, which requires 
that 

b± 2 = -b y2 3r / /21/2 • (19) 

Similarly demanding continuity of higher derivatives of p(r) will determine b i+1 , 2 for % > 2. 

Once A(r) and B(r) are determined from the two p(r) equations, p(r) is determined from 
the third equation. The three curvature invariants are also determined. We find that these 
latter four quantities contain a \r— r*| 1 / 2 e(r— r*) term with coefficients ^by 2 x [ai/r*, — ai/r*, 
— 2ai/r^, afbi/r^} respectively, and a |r — r*\ l l 2 term with coefficients |^ 2 X [ a i/ r *' a i/ r *) 
2ai/r^, — af^i/r*]. Given that the aj's and 6j's are determined from the two p{r) equations, 



5 This is more clear if one changes the coordinate from r to the proper distance from the crossing point 
oc \r — r*! 1 / 2 . 
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we conclude that there are two parameters by 2 remaining in the solution expanded around 
the crossing point. If these parameters both vanished then we would revert to the completely 
regular solutions in (JT3J). 

We now turn to the 2 + 4 gravity theory. A regular p(r) in this case does not imply a 
completely regular 1/A(r). To satisfy the various criteria above it turns out that we can set 
by 2 = a\ = = 0. Proceeding as before we find that continuity of p(r) fixes a 2 and b 2 
in terms of a^ 2 and b^, 2 . This then ensures that OR is also continuous. Continuity of p'{r) 
then determines a^, 2 and b^ 2 in terms of b^ 2 and a 3 , 2 . The coefficients of higher fractional 
powers are determined similarly. Thus in 2 + 4 gravity, the solutions of interest are described 
by the six parameters ai, £>i, b^ 2 and aj, 2 . If the latter four vanish, then we revert to the 
completely regular solutions. 

When we compare the solutions in GR and 2 + 4 gravity, satisfying the criteria above, we 
see that in 2 + 4 gravity the leading fractional power in the metric is 3/2 rather than 1/2. 
The same is true for p(r), so only in 2 + 4 gravity is p'(r) finite and continuous. Similarly 
OR is not continuous in GR, but it is in 2 + 4 gravity. Thus the transition to Euclidean 
signature has become smoother. In GR we have noticed that the relative size of |r — r*! 1 / 2 
and |r — r*| 1 ^ 2 e(r — r*) terms in p(r) and in the curvature invariants are determined by 
by 2 /by 2 . In 2 + 4 gravity each of the three curvature invariants can depend on a different 
combination of the same two terms, which depend in more complicated way on the four 
quantities by 2 and af/ 2 - In 2 + 4 + 6 derivative gravity we would expect that solutions 
satisfying the criteria above have b^ 2 = af, 2 = 0, and thus the curvature invariants would 
have a leading fractional power of 3/2 rather than 1/2. 

Even though the curvature invariants remain finite as the crossing point is approached, 
the existence of fractional powers in the metric can imply 5-function singularities in the 
curvature invariants. This can be checked independent of the field equations, which we have 
explained are ambiguous at r*. The three curvature invariants have (^-function singularities 
when by 2 are nonvanishing. These singularities are thus absent in our solution for 2 + 4 
gravity. But in this theory a 5-function singularity would arise in OR for example, due to 
the presence of b 3 , 2 and a 3 , 2 . In 2 + 4 + 6 gravity these would also be absent, and so on. 

Thus the basic pattern is that the leading singularities in higher derivative theories occur 
in higher derivative quantities, and in this way higher derivative theories become more and 
more smooth. This is the case as long as the higher derivative theory has a solution in which 
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more of the coefficients in ()17j) vanish. We must now confirm that such solutions exist, by 
showing that the solutions we have been describing near r* can actually be matched onto 
solutions that are regular at the origin and the horizon to form a complete black hole solution 
(the third criteria above). For this we need a numerical treatment. We also need a numerical 
treatment to check that the weak energy condition can be satisfied in the 2 + 4 theory. And 
we have already mentioned that it was the numerical study that prompted us to consider 
seriously the type 2 solutions, due to the difficulty of realizing type 1 solutions. 

VI. NUMERICAL STUDY 

For the numerical study the regular solutions at the origin and at the horizon with p(r/ l ) = 
v{ r h) = provide boundary conditions. But the equations cannot be numerically integrated 
starting right at the horizon. One strategy is to have the interior matter distribution vanish 
just before the horizon is reached, so that the Schd solution can be used to supply the initial 
conditions at a point just inside the horizon. The matter distribution would be continuous 
but nonsmooth at that point. Alternatively one can use a smooth matter distribution that 
vanishes at the horizon and has a only a small and high order leading derivative there. Then 
it can be an arbitrarily good approximation to use the Schd solution to supply the initial 
condition at a point just inside the horizon. There is a similar numerical problem with 
integrating from r = in 2 + 4 gravity. Here we use the series solution in (J2J) to provide 
initial conditions at a small but finite value of r. The solution then depends on the free 
parameters in this expansion, of which there are two in 2 + 4 gravity. 

For GR there are no undetermined parameters in the boundary conditions, once is 
given. In fact one of the equations can be solved analytically for 1/A(r) given p(r), and the 
constraints on p(r) can be satisfied even before B(r) is found. The numerical results for 
B(r) then confirms what we have already alluded to for GR: b^ 2 are in general nonzero. 
This eliminates the possibility of having the completely regular solutions in (fPfj) in GR, at 
least for the form of we are considering. 

Two examples of solutions are displayed in Figures (j2J) and 0- It is clear that we have 
numerically realized the criteria of the last section, where discontinuities are absent even 
though by 2 0. In particular b± is taking on the appropriate value, although the presence 
of this term is not directly visible in the B(r) curve. We see the presence of terms with 
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FIG. 2: A solution in GR with horizon size = 20, in Planck units. A) An infinite slope appears 
in B{r) at the crossing r* where 1/A(r) changes sign. B) The nonsmooth behavior of p(r) at the 
crossing is barely evident. C) A curvature invariant is compared to the same quantity for the Schd 
solution (in green). D) The region around the crossing point in (C) is enlarged. 

fractional power 1/2 in p(r) and the curvature invariants. These terms are much more 
evident in B(r) than in the physical quantities. The two sets of figures also illustrate how 
solutions exist for a wide range in the amount of matter present, and that the matter is 
more concentrated in the Euclidean core when more matter is present. 

We turn to the 2 + 4 theory. Here one must integrate the equations up from the origin and 
down from the horizon, and then manually adjust p(r) and the free parameters a.2 and 62 
appearing in 0, in order obtain solutions that match in an acceptable way at the crossing 
point r*. We verify the existence of solutions of type 2 for a range of the parameters a and 
(3 in the action (J3J. An example is displayed in Fig. (jlj) for a = —3/4 and (3 = 1. The 
results for A(r) and B(r) indicate the difficulties in numerically integrating the equations 
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FIG. 3: Another solution in GR with horizon size = 20, differing from Fig. © by having much 
smaller p and p. 

all the way to the crossing point, although physical quantities like p(r) are less sensitive to 
the numerical problems. We find that it is possible to realize the criteria of the previous 
section with b l > 2 = a\ = = 0. In particular we find that p(r) and its first derivative are 
both finite and continuous at the crossing point, so that its leading fractional power is 3/2 
rather than 1/2. Numerically we confirm that the curvature invariants are continuous, and 
in particular we find that for R the \r — r*| 1 / 2 e(r — r*) term dominates (see Fig. (jlP)) as it 
does for GR, while for R^R^ and R IMU(Tp R fJ ' UCTp the \r — r*! 1 / 2 dominates instead. 

In addition to finding that the transition to a Euclidean core is smoother in 2 + 4 gravity, 
from Fig. Q we see that the weak energy condition can be satisfied. We also see that 
the dominant energy condition is not satisfied in this case. At the origin of a Euclidean 
core we derived earlier the general constraint p(0) = — p(0), implying an approach to either 
Euclidean dS or AdS space near the origin. We find AdS for GR and dS for 2 + 4 gravity. In 
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FIG. 4: A solution in the 2 + 4 derivative theory of gravity with horizon size = 20, in Planck 
units. A) The metric functions (numerical problems are evident). B) The region around the 
crossing point is enlarged. C) p(r) remains positive and it and its first derivative are continuous, 
unlike GR. D) The scalar curvature is continuous and finite, as in GR. 

GR we find that the value of r* is quite flexible, while in 2 + 4 gravity we tend to find that 
r* remains small, not too much larger than the Planck length, even as gets large. 6 In this 
sense the solutions of the 2 + 4 theory are more restricted, although there still appears to 
be considerable freedom in the amount of matter in the black hole of a given radius. 

We have also done some exploration in 2 + 4 + 6 derivative theories. In these theories there 
are corrections to the Schd solution in the vacuum, and so the effects of matter are more 
difficult to disentangle. But we have confirmed the general tendency for a type 2 crossing 
point to form for some range of parameters in the action. There is even an indication that 
the crossing point can exist without the presence of matter in the region r* < r < ?+. 

6 This assumes that a and (3 are of order one. 
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When we compare Fig. (£Q) to Fig. (J2HH1) we see that the behaviors of p(r) and p(r) 
are distinctly different even in regions where the curvatures are small. To understand this 
we consider the deviation of the metric from the Schd solution in the presence of matter, 
expanded about the horizon. 

j X 

~rr~\ = l- — +J2ai(r-r h y 
Air) r 

K ' i=i 

oo 

B(r) = l-^+J^kr-r^f (20) 
i=i 

For GR b\ is a free parameter, for 2 + 4 gravity b\, i>2, Si are free parameters, etc. We can set 
these free parameters to zero to get as close as possible to the Schd solution. Then in 2 + 4 
gravity the deviation from the Schd solution in the presence of matter starts at higher order 
than it does in GR. In this case solutions to the 2 + 4 gravity equations close to the horizon 
do not even approximately satisfy Einstein's equations, even when the horizon radius is 
large. 7 It may seem odd that higher derivative terms, whose effects should be suppressed 
in regions of small curvature, are able to cause a breakdown of the Einstein equations. If 
we let 5A and SB denote the deviations from the Schd solution, the point is that at the 
horizon the derivatives of these quantities which appear in the Einstein equations are taken 
to vanish as initial conditions. Thus in the vicinity of the horizon the matter is influencing 
the higher derivatives of 5A and SB through the 4 derivative terms. For smaller r the 
curvature invariants gradually grow large and all terms and derivatives become of natural 
size. 



VII. DISCUSSION 



In conclusion we have searched for black hole solutions with a regular center. For a regular 
center to exist, at least one metric component must change sign in the interior of a black hole. 
Interior to this transition, spacetime may have Lorentzian or Euclidean signature. But we 
find that the Euclidean possibility is more typical and likely to occur. The requirements on 

7 We can consider the coefficients in (|2(J|) assuming that the leading nonzero derivative of p(r) at r = is 
p( n '(rh)- In GR we find that in the large limit that & n +i oc rhP^ n -'( r h) and b n — > 0. In 2 + 4 gravity we 
find that a n+ 2 — -~(2n + 3)b n+ 2 oc r^j^™' (rh) for n > 1. One result is that the value of Ricci scalar R in 
2 + 4 theory is larger than that in GR when 1/7-^ < \r — 1). 
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the energy momentum tensor at the crossing is less constraining, and the numerical analysis 
where one integrates the equations in from the horizon indicates that the Euclidean case is 
more generic. Although the weak energy condition is not satisfied near the crossing point 
in GR, it can be in higher derivative gravity. 

The possibility of a regular black hole raises the issue of whether such solutions can be 
continuously connected to regular solutions without a black hole, where A(r) and B(r) are 
everywhere positive. At the transition between these classes of solutions the crossing point 
would have to merge with the horizon, so that 1/A(r) vanishes at but remains positive 
elsewhere. The horizon then becomes extremal. For the black hole with Euclidean core B(r) 
is negative everywhere inside r^, and thus B(r) should vanish everywhere inside for the 
extremal solution. This allows it to continuously and smoothly connect to solutions where 
B(r) is everywhere positive. Interestingly just such extremal solutions have been presented 
for the Einstein- Yang-Mills theory in 11]. 

We have found that near the point of signature change the solution will generically depend 
on fractional powers of the distance from the crossing (in the coordinate system we have 
chosen). Nevertheless we have confirmed the existence of solutions where p(r) and p(r) and 
various curvature invariants are at least finite and continuous in GR. There are smoother 
solutions in the 2 + 4 derivative gravity theory where the leading fractional power is higher 
than in GR. We expect that the trend continues for theories with even more derivatives. 
Then as the number of derivatives in the theory increases, the curvature invariants with 
singular behavior will have to involve more derivatives. 

At the same time that the leading fractional powers in the metric increase, more of the 
terms that are regular (the ones appearing in ()14Jl ) are left undetermined by the equations 
expanded around the crossing point. The number of these undetermined parameters equals 
the number needed to match onto the physical Newtonian solution at large r. Thus we 
see how we are recovering the completely regular solutions in (j!4)l order by order in the 
expansion, as we consider theories with more and more derivatives. The implication is that 
a completely smooth transition to a Euclidean core emerges in a theory with all derivatives, 
ie. the complete derivative expansion of the true underlying theory. 

This picture is prompting us to contemplate the underlying theory. We are also prompted 
by the way our description patches together Lorentzian and Euclidean theories of classical 
gravity, leaving the description of the crossing point itself ambiguous. A theory that describes 
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a smooth d yna mical change of signature may require more degrees of freedom or higher 



dimensions 



12| . The notion that spacetime signature is not fundamental also emerges in 



string theory, where it is found that a string (M) theory on a Lorentzian torus is T-dual to 
a theory with a different signature . And finally, in the spirit of quantum cosmology, the 
Euclidean core can perhaps be viewed as a Euclidean instanton, representing a tunnelling 
amplitude from the black hole interior to nothing. 

In section 4 we described a smooth time dependent orbifold construction that could realize 
a regular black hole without a center. Here we should consider the more general solutions 
in section 5 where, for example in GR, we have 61/2 7^ 0. In the orbifold construction this 
implies that the metric depends on \u\, which in turn induces a (5-function in T^. Since there 
is no signature change, this 5-function now has some meaning. In other words, some kind 
of brane-like source is needed at the space-like boundary. In theories with more derivatives, 
where the leading fractional power is higher, it appears that a 5-function would arise in 
the highest derivative terms in the field equations, again implying a 5-function in T^ v . But 
as before, a singularity in any given curvature invariant would be absent in a theory with 
enough derivatives, meaning that the orbifold geometry is becoming smoother. Whether a 
completely smooth orbifold can be realized remains an open question. 

We have labelled some of our results as robust. In particular we have noted some prop- 
erties of the crossing point, the origin of the Euclidean core, and the event horizon which 
are true in GR, 2 + 4 derivative gravity, and 2 + 4 + 6 derivative gravity. We believe that 
these are examples of generic predictions of classical gravity. 
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